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Models of contagion arise broadly both in the biological and social sciences, with applications 
ranging from the transmission of infectious diseases to the diffusion of innovations and the spread 
of cultural fads. In this Letter, we introduce a general model of contagion which, by explicitly 
incorporating memory of past exposures to, for example, an infectious agent, rumor, or new product, 
includes the main features of existing contagion models and interpolates between them. We obtain 
exact solutions for a simple version of the model, finding that under general conditions only three 
classes of collective dynamics exist, two of which correspond to familiar epidemic threshold and 
critical mass dynamics, while the third is a distinct intermediate case. We find that for a given 
length of memory, the class into which a particular system falls is determined by two parameters, 
each of which ought to be measurable empirically. Our model suggests novel measures for assessing 
the susceptibility of a population to large contagion events, and also a possible strategy for inhibiting 
or facilitating them. 
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Defined broadly as the transmission of an influence 
from one individual to another, the concept of contagion 
occupies an important place both in biology — specifically 
in mathematical epidemiology 0, 0] — and in the social 
sciences, where it is in manifested in problems as diverse 
as the diffusion of innovations [E 0] , the spread of cul- 
tural fads |E IE j an d the outbreak of political Q or 
social unrest. 

Despite the wide range of social and biological phenom- 
ena to which they have been applied, existing models of 
contagion typically fall into one of two categories that we 
delineate in terms of the relationship between successive 
exposures of a "susceptible" to an "infectious" individual: 
(1) what we call "Poisson" models, in which successive 
contacts result in contagion with independent probabil- 
ity p; and (2) "threshold" models, in which the proba- 
bility of infection changes rapidly from low to high as 
a critical number of simultaneous exposures is exceeded 
(thus the effect of any single exposure depends strong- 
ly on the number of other exposures). The SIR mod- 
el 10], the canonical model of biological contagion, is an 
example of a Poisson model, as is the oft-cited Bass mod- 
el 3] from the diffusion of innovations literature. By con- 
trast, numerous models in sociology 0, economics |ll| . 
and political science 1121 . are explicitly threshold models; 
while others still |l3l Il4l llij embed thresholds implicit- 
ly through the relative costs and benefits of one action 
versus another. 

None of these models, however, treat the interdepen- 
dencies between exposures themselves as an object of 
study — rather they are simply assumed to either exist or 
not exist — hence their effects on the collective dynamics 
of contagion are unknown. Furthermore, if, as we show 



below, these effects turn out to be considerable, existing 
models provide no way to determine the conditions under 
which one kind of collective behavior or another should 
be expected. 

In this Letter, we explore a generalized model of con- 
tagion that, by introducing memory of past exposures 
to a contagious influence, generalizes and interpolates 
between Poisson and threshold models of contagion. Our 
model is defined as follows. Consider a population of 
N individuals, each of which is in one of three states S 
(susceptible) , I (infected) , or R (removed) . At each time 
step t, each individual i comes into contact with one oth- 
er individual j, drawn randomly from the population. If 
i is susceptible and j is infected then, with probability 
p, i receives a positive dose di(t), drawn randomly from 
some distribution of dose size f(d); otherwise, di{t) = 0. 
Each individual maintains a memory of doses received 
over the previous T time steps, recording a cumulative 
dose Di(t) = J2t'=t-T+idi(t')- Susceptible individuals 
become infected if Di(t) > d* where d* (i's dose thresh- 
old) is drawn randomly at t = from a distribution g(d*), 
and remains fixed thereafter. The probability that a sus- 
ceptible individual who encounters K < T infected indi- 
viduals in T time steps will themselves become infected 
is therefore 



fc=l v / 
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where 
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FIG. 1: Dose response curves for T = 10 (see Eq. Q): (a) Poisson (e.g., SIR-type) model: probability of receiving a positive 
dose by contacting an infective p = 0.3, distribution of dose sizes f(d) = S(d — 1), and distribution of individual dose thresholds 
g(d*) = S(d* — 1); (b) Stochastic threshold model: p = 1, f(d) is distributed lognormally with unit mean and variance 0.4333, 
and g(d*) — S(d* — 5); and (c) Deterministic threshold model: p = 1, /(d) = S(d — 1), and g(d*) = 5(d* — 5). 



is the average fraction of individuals infected after receiv- 
ing k positive doses in T time steps, and f(2i=i di > d*) 
is the probability that the sum of k doses drawn from f(d) 
exceeds a given d* . 

Equation Q can be thought of as an average dose- 
response relationship [lf3 | for the population in question. 
Figure ^ displays three examples of Eq. JQl for different 
choices of p, f{d), and g(d*). When all doses di — d are 
identical, all members of the population have the same 
threshold d* — d and p < 1, then Eq. JQl reduces to 
the standard SIR model, Fig. ^a); and when p = 1 and 
d* > d, it is equivalent either to a stochastic 0, fl5| . 
or deterministic 0, 0, threshold model, depending 
on whether doses are allowed to vary, Fig. ^b), or are 
identical, Fig. ^c). More complicated choices of f(d) 
and g(d*) correspond to many other kinds of models that 
incorporate varying degrees of interdependency between 
contagion events and also heterogeneity across individu- 
als. 

Once infected, individuals may recover with proba- 
bility r if Di(t) falls once more below d* (otherwise 
they remain infected), and recovered individuals become 
re-susceptible with probability p. While the resulting 
dynamics is, in general, quite complex, in the special 
case of p = 1 and r = 1 (analogous to so-called SIS- 
type [ija, LLi| models with instantaneous recovery), we 
can write down an equation for the steady-state fraction 
of infectives (fixed points) in the population |27| : 

T 

^ = E Qwfa-pfrft. (3) 

with Pfc as denned by Eq. J5J. 

Exact solutions of Eq. © may be obtained numerically 
(i.e., to arbitrar y p recision) and agree with results from 
our simulations 20]. Moreover, we are able to determine 
analytically that the equilibrium behavior of our gener- 
alized model falls into one of only three universal classes, 



examples of which are given in Fig. [21 We define these 
classes by the behavior of the fixed point curves around 
a single transcritical bifurcation |2lj ] , which is present for 
all models and located at p — p c — 1/(TP\) and <fi* = 0. 
We find that for a given value of the memory parameter 
T, the class into which a particular model falls depends 
only on the values of two quantities: Pi and P2, the prob- 
ability that an individual will become infected as a result 
of one and two exposures respectively. The three classes 
of behavior and their associated conditions are as follows. 

Class I: epidemic threshold models. When Pi > P2/2, 
we find the equilibria of Equation 2 fall on a bifurcation 
curve similar to that shown in Fig. |2fa), in which the 
sole bifurcation is the transcritical one at p = p c < 1 
(analogous to a continuous phase transition When 
P < Pc, a stable equilibrium exists at <p* — (i.e., all ini- 
tial seeds die out exponentially fast), and when p > p c , 
4>* > implying that a finite fraction of the population 
will become infected (i.e., an epidemic will occur). Thus 
the dynamics of models in class I is qualitatively equiv- 
alent to that of SIR-type models in which p c , sometimes 
called the epidemic threshold [2j| , determines the critical 
value of the infectiousness p required in order for an ini- 
tial seed of infectives to trigger an epidemic. We therefore 
call models in class I epidemic threshold models. 

Class II: vanishing critical mass models. When P2/2 > 
Pi > 1/T, we find, as shown in Fig. |2fb), a change in 
the nature of the transcritical bifurcation. Whereas for 
class I models, a stable fixed point curve emanates from 
(pa 0) with positive slope, class II models have an unsta- 
ble fixed point curve entering the transcritical bifurcation 
with a negative slope. Accompanying this is the appear- 
ance of a saddle-node bifurcation |21j at p = pb < p c 
and (f> = (f>l > (the system now exhibits a first-order 
phase transition). The solid (dashed) lines in Fig. Etb) 
correspond to the position of stable (unstable) equilib- 
ria respectively; thus there exists a region in p in which 
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FIG. 2: Example fixed point curves for the three basic classes of contagion models (all curves are obtained numerically using 
the closed form expression of Eq. © and agree with results from simulations (not shown) fx!^): a) Class I: Epidemic threshold, 
b) Class II: Vanishing critical mass, and c) Class III: Pure critical mass. Dose sizes are lognormally distributed with mean 1 
and variance 0.433, T = 10, and thresholds are uniformly set at a) d* = 0.5, b) d* = 1.6, and c) d* = 3. In the class II example 
of (b), trajectories of two initial conditions are indicated by arrows. For p = 0.445 and (f>(0) = 0.174, the contagion fails to 
persist while for the same p and a slightly higher initial level of infection, <f>(0) = 0.175, the contagion is sustained with 90% 
of individuals eventually infected. At the transition between classes I and II, the saddle-node and transcritical bifurcations 
coincide yielding a continuous phase transition distinct in nature from those of class I. The behavior of the fixed point curve 
near p — p c is (f>* cx (p — Pc) 1 ^ 2 rather than cf>* tx (p — Pc) 1 , though in highly special cases the exponent is lower |20||. Inset in 
(a): An example of a more complicated fixed point diagram. Here, T = 20, dose size is set to unity, f(d) = 5(d— 1), and d* = 1 
with probability 0.15 and 6 with probability 0.85. In principle, systems with more bifurcations are possible although unlikely 
since a strongly multimodal distribution of thresholds and/or dose sizes is required. 



two stable equilibria coexist, separated by an unstable 
equilibrium. In other words, as shown by the arrows in 
Fig. Efb) , if the initial seed <p(0) falls below the unsta- 
ble equilibrium, then the contagion will die out; whereas 
if it falls above the unstable equilibrium, then a large 
fraction of the population (corresponding to the upper 
stable equilibrium) will be infected (i.e., the system is 
metastable |22|). Thus infections in class II require a 
"critical mass" [24| to succeed. Because the size of this 
critical mass decreases to zero for p < 1, we call this 
class of models vanishing critical mass models, where we 
note that the sensitivity to initial conditions implicit in 
critical mass dynamics is absent entirely from models in 
class I; hence the two classes are qualitatively distinct in 
terms of their equilibrium behavior. 

Class III: pure critical mass models. When 1/T > Pi, 
the position of the transcritical bifurcation p c > 1; hence 
the only bifurcation potentially remaining in the interval 
< p < 1 is the saddle-node bifurcation, see Fig. |2fc). 
Thus in the manner of first-order phase transitions, a 
finite critical mass is always required (i.e., for all p) in 
order for an initial seed not to die out. As a result, we call 
models in this class pure critical mass models, where our 
classification of pure and vanishing critical mass models 
includes the dynamics of familiar threshold models 0,0, 
Il7j . but is more general. 

For some choices of f(d) and g(d*), additional equilib- 
ria are possible, and correspond to additional saddle-node 
bifurcations in Figure 2. Investigations of Eq. @ and 



numerical simulations |20j indicate that for additional 
equilibria to appear, one or both of f(d) and g(d*) must 
be multimodal, with widely separated peaks. For exam- 
ple, the inset of Fig.^a) shows the solutions of Eq. |J3J, 
where 15% of the population have d* = 1 and 85% have 
d* = 6. In addition to the standard Class I transcritical 
bifurcation, the inset shows a new saddle-node bifurca- 
tion. Extra saddle-node bifurcations generated in this 
manner, do not affect either the position or nature of 
the transcritical bifurcation, which continues to depend 
only on Pi and P2 . Hence our basic classification scheme, 
outlined above, is preserved for arbitrary f(d) and g(d*). 

The existence of three basic classes, and also sim- 
ple conditions that determine which class a particular 
contagion model belongs to, has interesting implications 
for understanding and possibly influencing the collec- 
tive dynamics of contagion, whether biological or social. 
First, the result dramatically reduces the effective com- 
plexity of the generalized model by showing that the qual- 
itative dynamics are independent of many of the details of 
the particular population [i.e., f(d) and g(d*)]; hence our 
use of the term "universal classes". Equally important, 
however, is that not all models fall into a single class; that 
is, not all kinds of contagion are the same. Furthermore, 
the three classes are non-degenerate in the rough sense 
that the conditions for each occupy broad regions of the 
parameter space. This result is not to suggest that some 
kinds of contagion are not more likely than others, but it 
does suggest that without empirical evidence about the 
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relevant Pi and P2, little can be concluded about the 
corresponding collective dynamics. 

Because neither the mathematical epidemiology liter- 
ature nor the micro-organismal dose-response literatures 
have tended to question the assumption that infection is 
a Poisson process, experimental evidence for or against 
the assumption is limited. In the context of biological 
contagion, therefore, our model suggests a very clear 
test to validate or refute the notion of interdependen- 
ces between exposures, and also specifies the strength of 
interdependencies required in order for them to be con- 
sidered important with respect to the collective dynam- 
ics. Specifically, if Pi > P2/2 then the assumption of 
independence between successive exposures can be con- 
sidered valid, whether or not it is precisely true, in the 
sense that the same qualitative behavior arises regard- 
less. If, however, Pi < P2/2 for some infectious agents or 
in some populations, then the interdependencies cannot 
be ignored, as their effect, in terms of the equilibrium 
states of the collective dynamics, may be dramatic. 

In the context of social contagion, our model has a 
different implication: that under very general, and one 
might argue likely, conditions, SIR- type models such as 
the Bass model that do not include memory or interde- 
pendencies between subsequent exposures, are incapable 
of capturing even the basic features of contagion dynam- 
ics. However, there are also likely to be some applications 
in which the independence assumption does turn out to 
be valid — that is, threshold models are also unlikely to 
be universally appropriate models of social contagion — 
where the key point is that our model provides precise 
conditions under which one state or the other will per- 
tain. 

A related point is that the conditions on Pi and P2 that 
separate the three classes of behavior suggest novel inter- 
vention strategies for suppressing, or alternatively stim- 
ulating, global contagion. Assuming that an individual's 
threshold d* can be manipulated — for example, increased 
(in a biological sense) by better preventative health treat- 
ment, or decreased (in a social sense) by exerting some 
financial, social, or cultural influence — then our results 
suggest that relatively minor manipulations (e.g., alter- 
ing Pi and P2 to shift p c while leaving the class of con- 
tagion unchanged) can have a dramatic impact on the 
ability of a small initial seed to trigger a global contagion 
event (i.e., the accessibility of a non-zero stable equilib- 
rium) as well as on the size of such an event if it occurs. 
We note that manipulating Pi and P2 is not equivalent to 
manipulating p since p determines the relevant point on 
the bifurcation curve, whereas Pi and P2 determine the 
shape of the curve itself; one cannot subsume temporal 
interdependencies within the independence assumption 
simply by lowering the per-event probability of transmis- 
sion. Finally, we note that while recent theoretical work 
on controlling disease epidemics, or by contrast stimulat- 
ing social contagion, has focused on individuals of excep- 



tional influence (e.g., so-called "super-spreaders" [25ll26j . 
and "opinion leaders" Ql), our model suggests that it 
could be the most easily influenced individuals (i.e., those 
contributing to Pi) who have the greatest impact on the 
dynamics of contagion. 
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